Dressed Wilson loops as dual condensates in response to magnetic and electric fields 
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We introduce dressed Wilson loops as a novel confinement observable. It consists of closed planar 
loops of arbitrary geometry but fixed area, and its expectation values decay with the latter. The 
construction of dressed Wilson loops is based on chiral condensates in response to magnetic and 
electric fields, thus linking different physical concepts. We present results for generalized conden- 
sates and dressed Wilson loops on dynamical lattice configurations and confirm the agreement with 
conventional Wilson loops in the limit of large probe mass. We comment on the renormalization of 
dressed Wilson loops. 



I. INTRODUCTION 

For the analysis of confining gauge theories Wilson and 
Polyakov loops are the key ingredients. Quark confine- 
ment reveals itself in an area law for the former and a 
vanishing expectation value for the latter. This is be- 
cause these loops represent non-Abelian factors along 
worldlines of heavy quarks. Consequently, the logarithms 
of the Wilson and Polyakov loop are related to the in- 
terquark potential and the free energy of a single quark 
at finite temperature, respectively. 

In a previous work, dressed Polyakov loops were intro- 
duced jT], which opened a new perspective on the rela- 
tion between confinement and chiral symmetry. Dressed 
Polyakov loops are center sensitive like conventional 
Polyakov loops and hence order parameters of the de- 
confinement transition [2]. Moreover, they connect to 
chiral symmetry since they are defined as Fourier compo- 
nents of quark condensates with respect to phase bound- 
ary conditions in the temporal direction reflected in 
the name 'dual condensate'. 

Dressed Polyakov loops contain a tunable probe mass 
m, which allows us to explore a number of their interest- 
ing features. In the limit of large mass one can expand 
the observable in inverse powers of the mass in terms of 
closed loops. Long loops are suppressed by the mass, 
which therefore prefers straight loops and makes contact 
with the conventional Polyakov loop in this limit. In 
turn, a finite mass governs the distance over which 'de- 
tours' in the dressed Polyakov loop extend. The probe 
mass also controls which part of the Dirac spectrum dom- 
inates the condensate. Furthermore, dressed Polyakov 
loops receive less renormalization than conventional ones 

a, cf. 0. 

In this work we introduce dressed Wilson loops defined 
as (collections of) planar closed loops of fixed area, but 
variable geometry. We view these loops as enclosing the 
same amount of disorder as conventional Wilson loops, 
which should reflect itself in an area law. In the same 
manner as for dressed Polyakov loops, dressed Wilson 
loops are defined using a tunable probe mass parameter. 
Large values of the probe mass can be shown to suppress 
long loops preferring 'ideal' loops, which on the lattice 
we demonstrate to be rectangular and thus make contact 
to conventional Wilson loops. 



Dressed Wilson loops are constructed via the quark 
condensate in the background of external abelian gauge 
fields. In this way, space-time and spatial Wilson loops 
are connected to the response of the quark condensate 
to (Euclidean) electric and magnetic fields, respectively, 
which gives an intuitive physical picture to be explored 
further. Similar magnetic fields play an important role 
at current heavy ion colliders and are investigated in re- 
lation to the QCD phase diagram. (Constant magnetic 
fields can be used to realize noncommutative coordinates, 
too, Wilson loops in this setting were calculated in |6j.) 

The definition of dressed Wilson loops through elec- 
tric and magnetic fields should also help to access these 
nonlocal objects in diagrammatic approaches and QCD 
models, as is the case for dressed Polyakov loops [?]• 

This paper is organized as follows. First we introduce 
external fields and define generalized and dual conden- 
sates/dressed Wilson loops, both in the continuum and 
on the lattice. 



After that, in Sect. Ill 



we present our 
numerical results from lattice simulations. Sect. IIVI is 
devoted to further conceptual remarks on our observ- 
able. In Sect. [V] we discuss the geometry of lattice loops, 
which enables us to make contact to conventional Wil- 
son loops in Sect. |VI| We conclude with remarks on the 



renormalization in Sect. [VTIl and an outlook. 



II. DRESSED WILSON LOOPS: IDEA, 
DEFINITION 

The concept of dressed Wilson loops is based on the 
possibility to incorporate the area of any (planar) loop 
in factors of the links along the loop. This is achieved by 
an auxiliary field strength, which is constant and abelian, 
through Stokes' theorem. 

To consider planar loops in two fixed directions /i and 
ly, let and a^, denote abelian gauge fields generating a 
constant field strength / with these indices, 

(!) 

ffj.i' — d^ajy — Syfl^ = / — constant for given /i, v (1) 

A possible choice on the infinite plane is a^, — fxp,,a^ — 0. 

For every given nonabelian configuration from the 
QCD path integral we modify the original nonabelian 
fields (indicated by capital letters) by 
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(in these directions only). Consequently, a closed loop 
on a contour C bounding an area S is multiplied by 



(3) 



which commutes through the nonabelian part (we set the 
charge to unity) . The area S is oriented (from the contour 
C), inverse loops are multiplied by the inverse factor. 
These modified links are to be used in a physical quantity, 
e.g. the condensate, and its dual quantity after Fourier 
transform in / will consist of Wilson loops of fixed area 
S. 

In order to specify the details, let us consider from 
now on the finite volume setting, i.e. the two directions 
being periodic with extensions and spanning a 



total area S 



Then the flux of abelian fields 



like momentum - obeys a quantization condition [5], 



fc e z. 



(4) 



in terms of the area, see e.g. [5]. 

We start with the generalized quark condensate 

S^j_y ^ Dk + m) 

where Dk is the Dirac operator evaluated with the ex- 
ternal abelian field of strength fk from Eqn. Q and 
restricted to the (/i, j/)-plane. The trace stands for a 
sum over color and spin degrees of freedom plus an inte- 
gral over space normalized by the two-volume S"^!, (such 
that tr is dimensionless and later on the lattice it simply 
amounts to a sum over all sites). Note that like for the 
four-dimensional condensate S is an intensive quantity, 
but with dimension of mass (not mass cubed). 

From gauge invariance it is clear that this observable 
consists of closed loops. Furthermore, the trace ensures 
that only connected loops contribute to Sfc (the loops 
are on the other hand allowed to self-intersect). Of these 
loops, we are interested only in non-winding loops (and 
only for those loops Stokes' theorem applies). We there- 
fore assume to have control over the contribution from 
winding loops, e.g. by a Fourier transform wrt. bound- 
ary conditions as for the dressed Polyakov loop (now in 
both compact directions). The role of winding loops will 
be elaborated on in more detail later. 

As argued above, closed loops in the generalized quark 
condensate obtain factors according to their area S G 
[0, S'^y] in the {fi, i/)-plane and Sfe is decomposed as 



dS t{S) e 



(6) 



This is nothing but a Fourier transform with the abelian 
fields fk and the area S as conjugate variables, which 
is in complete analogy to boundary condition angle and 
winding number for the dressed Polyakov loop. 

We invert this formula to obtain the dual condensate 
1 ^ 
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FIG. 1: The generalized quark condensate, Eqn. ([5|, in re- 
sponse to a constant magnetic field as a function of the flux 
quantum for several probe masses. On the upper axis we also 
show the magnetic field in physical units. 



which consists of loops having fixed area S. This defines 
the dressed Wilson loop. We constructed it to have di- 
mension mass cubed. The geometries of the contained 
loops are arbitrary, they contribute to Ti{S) with weights 
depending on m and the underlying space, see below. 

To calculate this novel observable on a discrete lattice, 
we have to specify the definitions a bit further. With 
lattice spacing a and A^^^^ — L^^^/a points in the two 
directions, respectively, the area of any loop is of course 



quantized as S 



sa 



with integer s € [0, N^^ = N^N^, 



The magnetic flux is also bounded from above. 



(8) 



(the value at which each plaquette receives a factor 1), 
such that all lattice quantities are A'^^-periodic in k. The 
external abelian fields can be implemented on the lattice 
by multiplying the link variables with appropriate space- 
dependent complex phases, see e.g. [TO]. 

The abelian factor for a Wilson loop of area s is 
exp{ifkS) = exp(27rifcs/A"^i,), where the flux and area 
quanta k and s are dual variables running over the same 
discrete range. The dual condensate/dressed Wilson loop 
is again defined via a Fourier transform, now the discrete 
one 



(9) 



fe=i 



(such that S(0) is the average of S^'s divided by 
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FIG. 2: The (absolute value of) dressed Wilson loops, 
Eqn. (|9|, dual to the condensates of Fig. [l] together with 
conventional Wilson loops (the latter measured for different 
rectangles with side lengths up to Ns/2). 



III. LATTICE RESULTS 

For our lattice studies we used the Symanzik improved 
gauge action and Nf — 2 + 1 flavors of stout smeared 
staggered fermions with physical masses. We generated 
5 configurations (which due to self-averaging effects of the 
Wilson loops sufficed for our purposes) on 16^ x 4 lattices 
at /3 = 3.31 (a w 0.4 fm, T « 120 MeV). Further details 
of the simulation setup can be found in e.g. |llj . The 
condensates, Eqn. (5|, were calculated from the spectra 
of Dk, cf. Eqn. ([lO); we computed all its eigenvalues at 
all possible magnetic fields with LAPACK. 

We present results for the staggered Dirac operator and 
the case of a constant magnetic field in the z direction, 
i.e. for Wilson loops in the {x, j/)-plane. The generalized 
quark condensate as a function of the magnetic field 
is shown in Fig. [l] The condensate at fc = 0, the con- 
ventional one without magnetic field, increases with the 
mass (for small masses) as is well-known. For nonzero k 
we expect that heavy particles are less affected by exter- 
nal fields and indeed the curve flattens out. The qual- 
itative behavior of is similar if one uses the 4d Dirac 
operator in the definition (not shown, for this check we 
used random estimators to determine the condensate). 

The dressed Wilson loops S corresponding to the con- 
densates shown in Fig. [T] are plotted as a function of area 
in Fig. [2] They reveal the expected decay in area similar 
to conventional Wilson loops. Our observable contains all 
kinds of loop geometries and therefore has an enhanced 
signal/noise ratio compared to conventional Wilson loops 
(by up to an order of magnitude for low masses). The 
decay rate as well as some modulations in its dependence 
on area, however, renders it different from conventional 
Wilson loops, to which we will come back below. 

The presented data are based on all closed loops, i.e. 
also winding loops. Since the number of the latter is only 



proportional to the circumference of the lattice (while the 
number of all the loops is proportional to the area) we 
expect winding loops to give negligible contributions to E 
for small probe masses. To confirm this we repeated the 
measurements under all combinations of periodic and an- 
tiperiodic boundary conditions in the two directions. Av- 
eraging these results (which can be considered the crud- 
est approximation to the zeroth Fourier component) ex- 
cludes loops with odd winding numbers. We observe that 
the results are within statistical errors unchanged. 



IV. FURTHER CONCEPTUAL REMARKS 

The applied constant abelian field /^^ is sensitive to 
loops in the (/i, i^)-plane. This means that electric and 
magnetic fields correspond to space-time and spatial Wil- 
son loops, respectively. While at zero temperature the 
QCD vacuum is 4d-symmetric, at finite temperature a 
difference is expected in the loops and correspondingly 
in the response to those fields. 

The abelian field is not sensitive to the '4d area' of 
loops extending in directions other than (/x, i/). Such de- 
tours would largely spoil the identification of loops of 
same area. Therefore, since we want keep contact to 
conventional planar Wilson loops of some area, we have 
restricted all generalized and dual observables, in partic- 
ular the Dirac operator, to some fixed (/i, z^)-plane. The 
calculation of its spectrum on the lattice is thus very 
cheap. This, however, does not mean that the system re- 
duces to a 2d gauge theory, since the gauge fields in the 
observable have interacted via the full action and hence 
embody 4d dynamics. Naturally, the calculation in a 
fixed plane is followed by an average over all such planes, 
just as for conventional Wilson loops. 

We emphasize that for the computation of dressed 
Wilson loops all possible abelian fields need to be in- 
tegrated over in the Fourier transform. (Typical mag- 
netic fields of phenomenological relevance for heavy ion 
colliders are physically very large, but correspond in cur- 
rent lattice simulations to the smallest possible quanta 
k = 0(1) — 0(10).) We also stress, that our setting 
is partially quenched in the sense that observables are 
measured on configurations without external field in the 
fermionic action (which should explain some of the dif- 
ferences to measurements in e.g. [TO])- 

An intuitive physical picture of our construction can 
be drawn from Wilson loops extending in space and 
time (say in the (a;, i)-plane), which represent a quark- 
antiquark pair. The corresponding dressed Wilson loops 
are obtained from an external electric field (Ex) instead 
of a magnetic one discussed so far. In the common picture 
of confinement the quark and the antiquark are bound 
together by nonabelian forces, that - if one was able to 
separate them to a large distance R - result in a potential 
V{R) growing linearly with R. Our (Euclidean) electric 
field acts on the quark and antiquark, as they carry oppo- 
site electric charges, and thus indeed attempts to vary the 



4 



distance R between them. Therefore it is clear that the 
electric field probes the nonabelian forces. Our formal- 
ism says that the latter are exhibited from the response 
of QCD quantities to all such electric fields (similar to 
inverse scattering). 

The QCD vacuum is symmetric under inverting the 
external field / — — /, i.e. k — )■ —k. The finite lat- 
tice resolution and the associated flux periodicity lead 
to the symmetry E/j = Eat^^^-^, obeyed on average, see 
Fig. [T] In order to obtain dressed Wilson loops at large 
areas, the high frequency fluctuations in Ej. need to be 
computed. Their calculation is limited by the flux quan- 
tization in finite volumes. The finite volume also limits 
large conventional Wilson loops. 

In the spectral representation of the Dirac operator 

(i) 

Dk with eigenvalues iAj, , that are purely imaginary and 
come in complex conjugate pairs (neglecting nongeneric 
zero modes). 



tr 



1 



D, 



E 



E 

j,A<''>0 



2m 



A 



W2 



(10) 



it is obvious that the dominant contribution to all con- 
densates comes from the IR part A < m (for the analo- 
gous discussion in dressed Polyakov loops see 12 ). 

An important effect in dressed Wilson loops is the sup- 
pression of long loops by heavy probes. For large m one 
can expand the condensate into a geometric series. 



tr ■ 



1 

m ^-^ 

1=0 



(11) 



For a lattice Dirac operator hopping over just nearest 
neighbors like the staggered operator we use, I is the 
length of the loop in lattice units, the length in physical 
units being L = la. To generate a closed loop, I needs to 
be even and the sign on the right hand side of Eqn. ( 11 ) 
can be removed. 



From Eqn. ( 11 ) it is clear that loops will be suppressed 



by the associated factor of inverse m (the same argument 
would apply to the scalar propagator from the lattice 
Laplacian). Thus we expect dressed Wilson loops to be 
'fuzzy' with the inverse mass governing their width. We 
will demonstrate this now on the lattice. 



GEOMETRY OF LATTICE LOOPS 



denotes the absolute value of the dual condensate (this 
also applies to the dual power of the Dirac operator T, 
defined later). 

What is the 'ideal lattice loop' in the sense of max- 
imizing the area at given circumference? For every 
concave loop the hooks can be folded to the outside 
keeping or even lowering the circumference. 




and the limits of this procedure are rectangles (among 
rectangles with same circumference those closest to 
squares have maximal area), contrary to the continuum, 
where the ideal loops are circles. 

The 'entropy', the multiplicity of different loops in the 
observable, may counteract the suppression away from 
the ideal loop. To check for this effect we define a com- 
binatorial factor 

, # loops of area s and circumference I , 
F[sJ) = (12) 

Note that the loops can start at any point in the lattice 
plane, therefore we divide by the lattice area to arrive at 
an intensive quantity. Furthermore, the circumference is 
meant in a general sense as the number of lattice bonds 
needed to traverse a loop of area s (since the link can, 
e.g., go back and forth). 

We calculated F using a recursive algorithm. This 
combinatorial factor apparently arises in different con- 
texts [12l Us]. At low area we reproduced the re- 
sults of |13| and extended them up to higher circumfer- 
ences/areas, see the appendix. At fixed area and large 
circumference F grows exponentially. 



Fis,l) 



for ; > s 



(13) 



The increase in this multiplicity factor (at fixed area s) is 
always exceeded by the increase in the suppressing mass 
factor m' for large enough mass (note that the color fac- 
tor from the loop depends on too, but cannot grow 
exponentially). Thus we conclude that the dual conden- 
sate E(s) at large mass contains predominantly loops of 
small circumference. 



In the following we specialize to the staggered Dirac 
operator again. The signs in this operator yield signs 
in the dressed Wilson loops depending on their length 
and area as follows. A closed loop consists of an even 
number / of links, of which 1/2 are inverse links W each 
coming with a minus sign. For the total staggered phase 
one can tile the area of the loop into plaquettes and con- 
vince oneself that odd areas get an additional minus sign, 
altogehter giving rise to a sign factor of (—1)'*+'/^. In or- 
der not to complicate the notation, in the equations E 



VI. CONTACT TO CONVENTIONAL WILSON 
LOOPS 

From the considerations in the last sections the dressed 
Wilson loops in the large mass limit should be propor- 
tional to conventional Wilson loops. Note that in this 
limit the UV eigenvalues now dominate, see also |14) . 

From Eqn. (11) it is clear that the leading contribu- 
tion in E(s) comes from the Dirac operator raised to 
the minimal circumference for the given area, denoted 
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by ?min(s), in other words, Zinin(s) is the first number I of 
link hoppings that is able to generate an area s (see also 
the appendix). Loops with areas s and circumference 
^min(s) describe certain geometries. For square number 
areas s — these are square loops p x p. The starting 
point can be anywhere on the square and therefore the 
combinatorial factor F is the number of points on the 
square 



F{p^4p)^4p 



(14) 



For areas s = p{p+l) the minimal loops are rectangles px 
(p + 1), where the longer side can point in both directions 
on the plane. Accordingly, the combinatorial factor is 
twice the number of points 



F{p{p+l),Ap + 2) = 8p + 4 



(15) 



In other cases the geometry of loops with minimal cir- 
cumference is not unique, but always 'close to rectangu- 
lar'. 

For powers of the Dirac operator we introduce the 
shorthand notation, 



T. 



1 



5, 



(16) 



which is an intensive quantity, and its Fourier transform 
T''(s) is obtained like in Eqn. Then for large probe 
masses 

where the correction term comes from the loop with same 
area but circumference increased by 2. This correction 
enters both the power of the Dirac operator and the mass 
suppression factor. 

For Dirac operators containing only U/2a and JJ*^^-* /2a 
one has the following relation to Wilson loops 



1 -F(s,Un(g)) 

a4 (2a)'-»('') 



N,{W{s)) (18) 



The color factor Nc comes about in order to take into ac- 
count the 1/Nc prefactor in the usual definition of traced 
Wilson loops W. 

Specializing to square loops, s = p^ with ^min = 4p and 
F from Eqn. ( 14 1 we obtain 



Combining that with Eqn. ( |17[ ) and reverting to physical 
circumference ('length') of Wilson loops L = Apa and 
side length R = pa = L/4, conventional square Wilson 
loops can be recovered from 



^E(i?2) ^ -1 ^ , NJW{RxR))(i+0( , ^ 

^ ' a4 (2am)^/'^ ^ ^ ^^V \{2am)^ 

(20) 
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FIG. 3: The leading term in the large mass limit of the dressed 
Wilson loop, (absolute value of) T'"'" (s), rescaled by the 
according combinatorial factor F and the other factors in 
Eqn. (18 1, in comparison to conventional Wilson loops. 



where the first factor on the right hand side accounts for 
the dimension of the dual condensate and the last one 
indicates the first correction from non-ideal loops being 
suppressed in powers of 2am (with coefficient that de- 
pends on the length). 

In order to prove the equivalence of the large mass 
limit of dressed loo ps a nd conventional loops numerically, 
j-'minCs) fj-oni Eqn. (18) is better suited than mS because 



the latter is tiny due to the according mass factor, see 
Eqn. (17). We show in Fig. |3] the numerical results for 
j-iminCsT^om Eqn. ( 18 ) divided by the accompanying fac- 



tors (note that the mass has dropped out of this relation) 
and see that indeed the conventional Wilson loops are re- 
produced. We note that one can prove that in this limit 
the contamination from winding loops becomes critical at 
s = {Ns/A)"^ + Ns/4: + 1 in the sense that winding loops 
contributing to this s can be shorter than ?min(s)- By 
filtering out loops with odd winding number once again, 
this critical area can be pushed up beyond / 4, where 
conventional Wilson loops are also out of reach. 



VII. REMARKS ON RENORMALIZATION 

From the connection of dual condensates to conven- 
tional Wilson loops at fixed lattice spacing it is tempting 
to obtain the renormalization of the latter^ from the for- 
mer, too. 

First of all, ordinary condensates defined with the 4d 
Dirac operator are known to be subject to additive and 
multiplicative renormalization. Dressed Wilson loops 
consist of subsets of loops present in 4d condensates, but 



^ The renormalization of Wilson loops is usually done via particles 
moving on them 1151 and is sensitive to edges of the Wilson loop. 



6 



feel the same lattice spacing and therefore should con- 
tain the same UV-divergences. For dual condensates one 
then expects that additive renormalization is not neces- 
sary since these divergences do not depend on the ex- 
ternal field^ and therefore are removed by the Fourier 
transform (as long as s ^ 0). The multiplicative renor- 
malization can easily be done by multiplying the conden- 
sate by the bare mass m, which has already been done 
in the previous section. 

Under the conjecture that mE(s) is a renormalized 
quantity, Eqn. (20) seems to reveal the a-dependcnt fac- 



tors needed to renormalize conventional Wilson loops'^. 
We remind the reader that this connection to Wilson 
loops is valid only for large am, which is of course related 
to the use of unphysical infinitely heavy probe quarks in 
conventional Wilson loops. In the continuum limit of 
condensates at fixed physical mass, however, the bare 
mass needs to run such that the combination am goes to 
zero PJ3j. More work is needed to understand this sit- 
uation (e.g. whether it also appears in other regulators 
than lattice) and the role of dressed Wilson loops. 



VIII. CONCLUSIONS AND OUTLOOK 

We have presented the definition and first results for 
dressed Wilson loops. They indeed obey a confining area 

^ We remind the reader that in our approach the external field 

only enters on the level of observables. 
^ and indeed the length dependence of these terms is of the form 

obtained in I15| . 



law and in the limit of large probe mass reproduce con- 
ventional Wilson loops. To demonstrate this we have 
used lattice simulations, but we strongly expect our ob- 
servable to be useful in continuum approaches as well. 

Being constructed from the response of quark conden- 
sates to magnetic and electric fields, dressed Wilson loops 
should help to get more insight into the interrelation of 
these concepts. This concerns in particular the behav- 
ior at the QCD transition. Further interesting aspects 
of dressed Wilson loops, partly deduced from the anal- 
ogy to dressed Polyakov loops, are their IR dominance 
and renormalization properties. Moreover, the possibil- 
ity to use dressed Wilson loops for setting the scale is 
attractive. 
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Appendix: Combinatorial factors 



In the following table we list the combinatorial factor F{s, I), the number of possible loops with I two-dimensional 
lattice hoppings ('circumference', but hoppings can be forth and back) that span an area of s lattice units, divided 
by the total number of lattice points, cf. Eqn. ( [T2| , up to I = 24. Note that the area is oriented and thus can be 
negative, the combinatorial factor is the same, Fj—s, I) = F{s, I). 



s\l 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 


22 


24 





4 


28 


232 


2156 


21944 


240280 


2787320 


33820044 


424925872 


5486681368 


72398776344 


972270849512 


1 




4 


72 


1008 


13160 


168780 


2168544 


28133728 


369612648 


4920045800 


66324542240 


904584355488 


2 






12 


308 


5540 


87192 


1291220 


18569808 


263462220 


3718483560 


52450578224 


741275064780 


3 








48 


1560 


33628 


610232 


10127744 


159762240 


2445203460 


36746804160 


546205048128 


4 








8 


420 


11964 


262612 


5015108 


88145244 


1469700900 


23667392012 


372270358704 


5 










80 


3636 


101976 


2289760 


45306288 


827935484 


14348361544 


239687474784 


6 










20 


1200 


40376 


1036368 


22761228 


452521560 


8408087996 


148916557340 


7 












264 


13720 


435040 


10920456 


239017700 


4784523128 


90018016224 


8 












72 


4900 


184104 


5208372 


124781340 


2683439616 


53545392516 


9 












12 


1512 


73056 


2398752 


63687440 


1479475008 


31384157616 


10 














420 


28064 


1085724 


32101920 


807111536 


18215426820 


11 














112 


10336 


478008 


15887160 


434604632 


10465152864 


12 














28 


3760 


208728 


7803420 


232384020 


5973330680 


13 
















1088 


84312 


3699300 


121904288 


3369519456 


14 
















352 


34560 


1749200 


63575468 


1889445720 


15 
















96 


13392 


805280 


32649760 


1048976880 


16 
















16 


4788 


359760 


16509284 


576770520 


17 


















1584 


155560 


8195968 


313483392 


18 


















540 


66960 


4037528 


169191804 


19 


















144 


26920 


1931952 


89921136 


20 


















36 


10680 


914540 


47383128 


21 




















3740 


415888 


24503088 


22 




















1360 


188144 


12567024 


23 




















440 


82016 


6330912 


24 




















120 


34496 


3138264 


25 




















20 


13640 


1517328 


26 






















5236 


723264 


27 






















1936 


337104 


28 






















660 


153360 


29 






















176 


66528 


30 






















44 


28464 


31 
























11280 


32 
























4488 


33 
























1632 


34 
























528 


35 
























144 


36 
























24 



The bold numbers are those with s = p or s = p{p + 1) and / = ^min(s) used in Eqns. (14) and (15). The geometries 
of the associated loops are ideal in that they maximize the area at given circumference (therefore in the table these 
numbers are the lower ends of the /-columns), namely by squares and rectangles, respectively. 
Eqn. ( 13 ) on the other hand describes the exponential growth of F in the rows at fixed area s. 

I 



